NASA-CR- 173445 
19840020197 


A Reproduced Copy 

OF 



***«»*«» no* 


this ROOJ4 


Reproduced for NASA 
by the 

NASA Scientific and Technical Information Facility 


LIBRARY COPY 

-■■VO 1984 


LANGLEY RESEARCH CENTER 
LIBRARY, NASA 
HAMETOH VIRGINIA 


i 


FFNo 672 Aug 65 


IQ 



OLD DOMINION UNIVERSITY RESEARCH FOUNDATION 




DEPARTMENT OF GEOPHYSICAL SCIENCES 
SCHOOL OF SCIENCES AND HEALTH PROFESSIONS 
0L0 DOMINION UNIVERSITY 
NORFOLK, VIRGINIA 


Technical Report No. GSTR-84-4 


IDENTIFICATION OF AEROSOL COMPOSITION 
FROM MULTI-WAVELENGTH LIDAR MEASUREMENTS 


By 


Sidney Allen Wood 


Technical Report 



Prepared for the 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, Virginia 23665 

Under 

Research Grant NCCI-28 4 

Dr. Edward V. Browell, Technical Monitor u 

ASD - Chemistry & Dynamics Branch 



( NAS A -CB- 1734 45) IDENTIFICATION OF AEEOSQ1 B84-28266 

COMPOSITION FEOfl HOLT I- H A V El ENGTH LIDAR 
NEASOBEBENTS (Old Doaiuiofl Dsiv., Mosfolk, 

Va.) 116 p HC A06/HF AOI CSCL 13B Unclas 

G3/45 00771 


April 1984 


5 * 


* 


- . /" 




f 


DEPARTMENT OF GEOPHYSICAL SCIENCES „ 
SCHOOL OF SCIENCES AND HEALTH PROFESSIONS 
OLD DOMINION UNIVERSITY 
NORFOLK, VIRGINIA 


Technical Report No. GSTR-84-4 


IDENTIFICATION OF AEROSOL COMPOSITION 
FROM MULTI-WAVELENGTH LIDAR MEASUREMENTS 


By 


Sidney Allen Wood 


Technical Report 


Prepared for the , . 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, Virginia 23665 


Under 

Research Grant NCCI-28 

Dr. Edward V. Browell, Technical Monitor 

ASD - Chemistry & Dynamics Branch 


Submitted by the 

Old Dominion University Research Foundation 

P.0. Box 6369 

Norfolk, Virginia 23508 


ABSTRACT 


IDENTIFICATION OF AEROSOL COMPOSITION 
FROM MULT I- WAVELENGTH LIDAR MEASUREMENTS 


Sidney Allen Wood 
Old Dominion University, 1983 
Director: Dr, Scott T. Shipley 


This paper seeks to develop the potential of lidar 
for the identification of the chemical composition of at- 
mospheric aerosols. Available numerical computations 
suggest that aerosols can be identified by the wavelength 
dependence of aerosol optical properties. Since lidar can 
derive the volume backscatter coefficient as a function of 
wavelength, a multi-wavelength lidar system may be able to 
provide valuable information on the composition of aero- 
sols. This research theoretically investigates the volume 
backscatter coefficients for the aerosol classes, sea-salts, 
and sulfates, as a function of wavelength. The results 
show that these aerosol compositions can be characterized 
and identified by their backscatter wavelength dependence. 

A method to utilize multi-wavelength lidar measurements to 
discriminate between compos itionally different thin aerosol 
layers is discussed. 
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INTRODUCTION 


Lidar (laser radar) techniques have been developed 
over the last two decades for a wide range of studies of 
atmospheric properties (Russell, et al., 1981). This paper 
seeks to broaden these lidar capabilities to include the 
identification of the chemical composition of atmospheric 
aerosols. Available numerical computations suggest that 
broad classes of aerosols can be identified by the wave- 
length dependence of aerosol extinction and backscatter. 
Since backscatter signals can be obtained by lidar as a 
function of wavelength, a multi-wavelength lidar system 
may provide valuable information on the composition of at- 
mospheric aerosols. 

This thesis theoretically examines the wavelength 
dependence of volume backscatter coefficients for the 
tropospheric aerosol classes, sea-salts (NaCl) and sul- 
fates ( ( NH^ ) 2 ^ 0 ^ ) . Such a study must account for wave- 
length variations due to the natural variability of 
particle si 2 e distribution, refractive index, and relative 
humidity. The model analysis characterizes sea-salt and 
sulfate aerosol classes by demonstrating that the compo- 
sitions are uniquely related to their backscatter spectral 
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variations. The results support the idea of using compo- 
sitional differences in a multi-wavelength lidar system 
to identify aerosol layers by their backscatter wavelength 
dependence. These results also provide useful information 
on the optimum lidar wavelengths for discrimination between 
the aerosol types considered. 

Previous model calculations and experimental lidar 
measurements suggest that the wavelength dependence of the 
aerosol volume backscatter coefficient varies with aerosol 
chemical composition. The capabilities of multi-wavelength 
lidar for the identification of various aerosol composition 
from characteristic backscatter signatures has been ex- 
plored, with moderate success, at infra-red wavelengths. 

The characteristic visible backscatter wavelength depend- 
ence expected for the aerosols in this study were estab- 
lished from previous pre-existing Mie model calculations. 

Aerosol backscatter coefficients have been theoreti- 
cally examined in the infra-red wavelength region 7.0 to 
1*1.0 pm (Wright, 1980). Mie model calculations for 
basalt, H 2 S0 4 ( 9 ^%) , H^SO^ (50%), ammonium sulfate 
(NH 4 ) 2 S0 4 , water and ice indicated significantly different 
backscatter features for each composition. Wright demon- 
strated that the six aerosol compositions were uniquely 
related by their infra-red backscatter spectral signa- 
tures. The author was also able to experimentally 
distinguish between high altitude ice clouds and low 



altitude water clouds with a multi-wavelength infra-red 
lidar system. 

Mudd et al. (1982) investigated the potential of 
using a multi-wavelength lidar to differentiate between 
two sulfuric species in the stratosphere. Volume back- 
scatter coefficients from ammonium sulfate and sulfuric 
acid particles were measured by a CO^ laser under labora- 
tory conditions, for the wavelength region 9.2 to 10.8 fim. 
Mie model calculations were made for comparison with the 
lidar measurements. The discrepancy between their experi- 
mental and theoretical values of volume backscatter coeffi 
cient wavelength dependence was suggested to be due to 
particle non-sphericity , inaccurate size distribution 
measurements or the presence of another species in the 
aerosol chamber. The author found that sulfuric acid and 
ammonium sulfate was distinguishable by backscatter signi- 
tures at infra-red wavelengths under laboratory conditions 

Available model calculations suggest that the broad 
aerosol classes, sea-salts and sulfates, can be identified 
by their backscatter wavelength dependence at visible 
wavelengths. Values of volume backscatter coefficients 
derived from the Mie calculations of Shettle and Fenn 
(1979), Dermendjian (1969), and Russell et al. (1981), in- 
dicate a wavelength dependence ranging from £ 2.0+0.32 
for sulfates to £ -0. 5+0.1 f or sea-salts. These calcula- 
tions were made for dry particles with one particle size 



distribution in the wavelength range 0,2 to 2.0 pm. 

This study describes two tropospheric aerosol models 
and their optical properties. The models were based on a 
literature review of available data on aerosol physical 
properties such as particle size distribution (aerosol par- 
ticle radius ) , complex refractive index (amount of scatter- 
ing and absorption) , and relative humidity. The optical 
properties of the aerosols were computed for the wavelength 
range 0.2 to 2.0 at various relative humidities ranging 
from 0 to 80 percent. The development of the two aerosol • 
models and their physical properties is presented in 
Chapter II. The computation of backscatter coefficients 
from aerosol optical properties at the small size para- 
meters 0 to 100 is discussed in Chapter III. Also included 
are the model calculations of backscatter wavelength de- 
pendence made for both sea-salts and sulfates. The dis- 
crimination between sea-salts and sulfates by backscatter 
wavelength dependence at various relative humidities is 
.presented in Chapter IV. A method to retrieve the wave- 
length dependence of volume backscatter from a two wave- 
length lidar system is discussed in Chapter V. 



II. MODELS OF AEROSOL PHYSICAL PROPERTIES 


Atmospheric aerosol models have been developed to in- 
vestigate optical properties of the maritime and continental 
classes of aerosols as a function of wavelength. Physical 
properties (i.e., particle size distribution, refractive 
index, relative humidity) for both of these aerosol types 
can vary significantly in the atmosphere, consequently in- 
fluencing the particle's optical properties. The aerosol 
physical properties used in the models are based upon mea- 
surements taken from the literature for various geographical 
locations, and seasons of the year. The measured range of 
natural variations found in both aerosol size distributions 
were incorporated into these models. Also included were the 
effects of increasing and decreasing relative humidity varia- 
tions on aerosol particle growth and refractive indices. 

The models for sulfate and sea-salt aerosols are presented 
in a way which attempts to account for natural variations in 
the particle size distribution, refractive index, and rela- 
tive humidity. 

2.1 MODEL SIZE DISTRIBUTIONS 

The aerosol size distributions used in the continental 
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and maritime models were based upon available experimental 
measurements, to provide a range of natural variability. 
Sulfate and sea-salt size distributions were collected from 
the literature for both airborne and surface locations, 
under various meteorological conditions. Some relevant in- 
formation (i.e., location, collection method and efficiency, 
meteorology, time of year) on the size distributions is pre- 
sented in Table la for sea-salts, and Table lb for sulfates. 
A normalization of the distribution curves revealed a sig- 
nificant size distribution variation in the small and large 
particle radius sizes as shown in Figs, la and lb. This 
natural variability of the particle size distributions was 
incorporated into the two aerosol models. 

A non-linear regression routine was used to best fit 
each empirically derived particle size distribution to a 
log normal curve. 


where 


dN(r) _ / l 
cJr vln (1 o) r<rv5tf, 


dN (r) 
dr 


)E XP |- (lo9r 2 4 0 ^ i 


is the total concentration of aerosol 
particles with radii greater than r, 

-3 

is the number of particles (cm ) , 
with radii between and r+dr, 


is the mode radius, 
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Table la. continued 
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2S><». for radiue O.Ol^a, the col- 
lection efficiency was 70%. 
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Table lb. Data on collected continental 
particle size distributions. 


Sfr e Cot Ik t Ion fethod 

Distribution IgciUw teteorploor end fffuieocy 

1 Sotittrfejrj foot, «t ry continental Aoyco tingle particle 

polar air eats. Starts counter. Collection 

♦atterljr winds. efficiency wi satis- 

Relative hualdlty ms factory for 0.4 * r< 

OS- 90S. Kaa. 

Z Socsterfeerg Cool, try continental Aoyco tingle particle 

polar air boss. »orth- counter. Collection 

easterly winds. efficiency was Mill* 

Relative hualdtty ns factory for 0.4* r< 

8S-90S. Uf*n. 

3 ffeppen Old earltlae air nets Royco tingle particle 

codified to be strongly counter. Collection 
continental. South- efficiency test satls- 

nsterty winds . factory for 0.4 1 r< 

Relative tewldlty tat l4/*a. 

SOS. 

4 Meppen Old nrfttne air asst Royco tingle particle 

sodHIed to t» strongly counter. Collect lor 
continental. South- efficiency m\ satis* 

westerly winds. factory for Q.A<r< 

Relative hwaidlty teat l4/»a. 

SOS. 

5 Roolos 1 leather was not given CxechotlavaMan neatrane 

Telep since the distribution filter 1.7a above growtd 

was average d over the facing wind direction, 

season. Counted and steed par- 

ticles In the radius 
range 0.0? -Q.S*«. A 
p&ase contrast optical 
•Icroscop? was used for 
particles O.S-lOO^r. 

The collection eff i 
cfency was 60S for 
0 . 0 S/<« 



(sulfate) 


ti»/ 

Al titude Source 

AM 17. If?* Cress, 1990 
1228 GSfT 
Alt: 1S0C * 


APR 12, 1978 Crest. 1990 

1248 car 

Alt: 5CCG» 


AUC S, 1977 Cress, 1980 
Alt: toils 


o o 

-T\ g 

-o o 

o r* 

o y 

r— 


<C 



r- 


2 i 


AUC S, 1977 Crest. I960 
947 cm 
Alt: 3060 b 


Suanrr fet zeros. 1|79 

1974-1978 

Alt: IXtOm 


t 



Table lb. continued. 


• Koalotl Heather mss not given Crechoslavafctan ««a*rene 

Telep since the <fitrtfc«tton filter 1.7w above jrnnd 

ms averted over the facing wind direction, 
season. Counted and sired par- 

ticles In the radius 
range 0.G7-0. *>M •- * 

phase contrast optical 
aicrosccpe was used for 
particles 0.5-100^n. 

The collection effi- 
ciency was 60s for 
0.0* ■. 


7 lost os 1 Heather ms not given CrechoslevaUan wewbrane 

Telep since the distribution filter t. 7a above ground 

MS averaged over the facing wind direction. . 

season. Counted and sired parti- 

cles in the radius range 
O.Q?-C.S^a. X phase 
contrast optical atcro- 
stGpe was used for par- 
ticles 0.S- 100 /»•- The 
collection efficient* was 
102 for r>0.0*^a. 

I 8/A* S/A tog norms 1 Curve fit to 

, KMtbr and Cantrell 

; (1976). Typical range 

of values for the 

li accumulation aerosol node 
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Fig. lb. Eight continental particle 

size distributions collected 
from available literature. 
The distributions vere 
normalized to 4 cm~3 for 
r=0.1/tm to show relative 
variability for the particle 
size distributions . 
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is the mode standard deviation 


The use of a log normal curve to represent particle size dis- 
tributions has been discussed in the literature (Whitby et 
al. , 1972; Kelkar and Joshi, 1975; Thomalla and Quenzel, 1982; 
Patterson et al., 1980). A comparison between a particle 
size distribution taken from Junge (1969) and the log normal 
curve used in the model to approximate the distribution is 
presented in Fig. 2. The general slope of the Junge dis- 
tribution is shown to be adequately simulated by a log normal 
curve even though there are some discrepancies between the 
two curves, especially at the larger radii. The values of 
the mode radius and mode standard deviation that best fit 
the distributions are shown in Table 2a for sea-salts, and 
Table 2b for sulfates. 


Table 

2a. 

Sea-salt 

Table 

2b. Sulfates 

SD 

r . 

l 

or 

l 

SD 

r . 

l 

<rr 

i 

1 

0.06 

0.40 

1 

0.07 

0.31 

2 

0.10 

0.37 

2 

0.05 

0.35 

3 

0.10 

0.33 

3 

0.05 

0.36 

4 

0.28 

0.21 

4 

0.07 

0.38 

5 

0.24 

0.25 

5 

0.07 

0.17 

6 

0.32 

0.39 

6 

0.09 

0.21 

7 

0.26 

0.44 

7 

0.07 

0.26 

8 

0.16 

0.40 

8 

0.07 

0.26 

F. 

X 

0.190 


r\ 

r 

0.0675 


n 

0.360 



0.2875 


°*i 

0.093 


°~x. 

1 

0.0128 


^ or. 

1 

0.092 


°*i 

0 ,0750 
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Fig. 2. A log normal size distribution 
best fit to a Junge, 1969 
particle size distribution. 
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The means and standard deviations of the log normal 
mode parameters (also shown in Tables 2a and 2b) were used 
with a random data generator to produce synthetic particle 
size distributions* The parameters r^ and <r^ were assumed 
to be normally distributed with variances as shown in 
Tables 2a and 2b. The randomly generated parameters which 
define the particle size distributions are presented in 
Table 3. 



Table 3. 

Mode Parameters 

for Synthetic Particle 

Size 



Distributions. 





5I7E 

MAR I T | mf 

(OCEANIC) 

CONTINENTAL (RURAL 


DISTRIBUTION r < 

o 1 

r 1 

° 1 


1 

. 355 

.470 

. OR 35 

. 374 


2 

. 163 

. ?95 

.0863 

. 231 


3 

. 340 

.363 

.0838 

. 303 


4 

.11? 

. 35? 

. 0679 

.278 


5 

.146 

. 364 

. 0638 

. 28 7 


6 

.?1K 

.44? 

. 07?7 

.352 


7 

. 190 

. 351 

. 0697 

.2 77 


8 

A 

. 1 70 

. 376 

. 0674 

.2 98 


M 

. 1 75 

.4 13 

. 06 79 

. 327 


1 0 

.16? 

. 366 

. 066? 

. 2«9 


1 1 

.195 

.310 

. 07R3 

. 243 


1 ? 

.147 

. 406 

. 064 A 

. 323 


1 3 

. 159 

. 3?8 

. 06R9 

. 258 

r 

1 4 

. ?46 

. ?53 

.0754 

.197 

r 

f ' 

1 S 

.114 

. 3 96 

.0583 

.315 

9 

I fi 

.149 

. 505 

. 0644 

. 403 


1 7 

.340 

. ?R4 

.0813 

. 222 


1 8 

. ?01 

.337 

.0710 

. 2*6 


1 9 

. ?50 

. 3P5 

. 0768 

. 305 


2 0 

. 1 96 

. 3 a 3 

. 07T4 

.27 1 

/ 

21 

.173 

. 353 

. 0676 

.279 


22 

. ? 5 1 

. 394 

. 0758 

.31? 


23 

. 156 

. 304 

. 0663 

.141 

* 

24 

.175 

. 36? 

. 0654 

. 238 

*f* 

25 

. ?70 

. 359 

. 0679 
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Since this study is mainly concerned with free tropo- 
spheric aerosol particles, only one mode (accumulation mode) 
was used for each distribution, even though size distribu- 
tions of aerosol particles have been documented as having a 
multi-mode structure (Whitby et al., 1972; Willeke and 
Whitby, 1975; Willeke, 1974). The use of only one mode for 
the particle size distribution at the higher altitudes is 
consistent with the results of Whitby and Cantrell (1975). 

2.2 THE EFFECTS OF INCREASING AND DECREASING RELATIVE 
HUMIDITY ON AEROSOL PROPERTIES 

In this section, the effects of humidity variation on 
aerosol properties and how they are incorporated into the 
models is discussed. As relative humidity increases, aero- 
sol particles in the atmosphere tend to grow as water vapor 
condenses onto the particle by the process of adsorption 
(Fitzgerald, 1975). The process will continue until at some 
critical relative humidity, the particle will undergo a 
sudden change to a saturated solution droplet by the pro- 
cess of absorption (Fitzgerald, 1975; Covert, Charlson, and 
Ahlquist, 1972). The process of particle water uptake not 
only increases the size of the particle, but also changes 
the aerosol chemical composition, and therefore the refrac- 
tive index. 

Numerous authors have addressed the problem of at- 
mospheric particle growth as a function of relative 
humidity (Orr et al., 


1958; Winkler and Junge, 1971 and 
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1972; Meszaros, 1971; Herbert, 1975; Fitzgerald, 1975; Tang, 
1976; Tanq et al., 1977a, 1977b and 1978; Thudium, 1978; 
Fitzgerald et al. , 1982) , and the foremost of the studies 
is the work of Hanel (1968, 1970, 1971, 1972a, 1976, 1979, 
1980) . Hanel constructed equilibrium growth curves of 
"typical" aerosol particles based on the characteristic be- 
havior of aerosol group samples. 

A semi-empirical model has been developed by Hanel to 
describe the growth of a particle's radius due to relative 
humidity. In order to derive the equilibrium growth of a 
particle, the water activity (i.e., the equilibrium relative 
humidity over a plane surface) of the solution must be mea- 
sured (Hanel, 1976). Hanel made the approximation that the 

volume (V) of a solution drop was equal to the size of the 

. 4 3 

volume of the dry spherical particle V 0 - fir,, , plus the 

A. 3 

volume of pure water V/j* -j- JY P, condensed on it 


V _ 

V.” 


7-f 


s UtJ Of 


( 2 ) 


The equivalent radius of a particle, the radius of a sphere 
with the same volume as the solution drop, can then be de- 
fined as 


=r,(i + 



(3) 
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where 


r 0 

P 0 


P w 


K, 


w 


M, 


is the equivalent radius, 

is the equivalent radius at a relative humidity 
of zero 

is the mean bulk density of the particle at 
RH=0 , 

is the density of water, 

is the mass of water condensed on the particle, 
is the mass of the particle in its dry state. 


The particle equivalent radius is reliable for calculating 
the humidity dependence on extinction and scattering coeffi- 
cients even at low relative humidities (Hanel , 1976). 

According to the refractive index measurements of 
Volz (1972, 1973), the amount of absorption of continental 
and oceanic type aerosols is quite low when compared with 
the amount of scattering at the visible wavelengths, there- 
fore only the real part of the complex refractive index was 
used in this study. The effective refractive index, m, can 
be easily calculated once the "wet" particle radius has 
been found by the mass weighted function by the equation 



(Hanel, 1976) 


(4) 
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where 

m is the real component of the complex refractive 
index for water, 

m is the real component of the complex refractive 
° index for the dry substance. 

The refractive indices for water as a function of wavelength 
were taken from the measurements of Hale and Querry (1973) , 
and are shown in Table 4. Also shown are the dry refractive 
indices for the maritime and continental models taken from 
Shettle and Fenn (1979) , which are based on the measure- 
ments of Volz (1972, 1972, 1973). 


Table 4. Refractive index of water, 
sea-salt, and sulfate. 
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The use of Hanel* s equations to describe the relative 
humidity effects on "typical” aerosol particles requires 
measurements of the mass ratio between dry substances and 
water, relative humidity, and the dry particle density. 

Such measurements have been made for maritime, industrial, 
continental, and desert type aerosols by Hanel (1976) and 
Hanel and Lehmann (1980) . Their data was collected for 
aerosol particles exposed to increasing and decreasing rela- 
tive humidity conditions. For decreasing relative humidity, 
a deliquescent material will contain more water than the 
particle would under increasing relative humidity. This re- 
sult is called the hysteresis effect and is mainly due to 
the phenomenon that evaporating salt solutions tend to 
supersaturate with respect to the salt in the solution. The 
hysteresis curves for maritime type aerosols have been found 
to be broader at higher relative humidities than continental 
type aerosols (Hanel and Lehmann, 1980), as shown in Fig. 3. 
These results suggest that maritime aerosol type particles 
can return to their original size at lower relative humidi- 
ties, whereas the continental type particles do not. 

2.3 MARITIME (OCEANIC) AEROSOL MODEL 

The composition for the maritime aerosol model is 


* 


& 


assumed to be pure sea-salt. Aerosol compositions through- 
out the marine boundary layer (MBL) and marine free 
troposphere (MFT) have been documented as being composed 



Fig, 3. Hysteresis curves for maritime aerosols 
(Mace, Ireland) and continental aerosols 
(Hohenpeissenberg , Germany) (Hanel and 
Lehmann, 1980). 
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of a mixture of oceanic aerosols and continental background 
aerosols (Junge, 1969; Meszaros and Vissy, 1974; Patterson 
et al., 1980). According to Meszaros and Vissy, the larger 
particles consist mainly of sea-salt, while the smaller par- 
ticles are the rural background aerosols. Particle size 
distribution measurements (Patterson et al., 1980) show that 
the larger sea-salt particles in the MBL predominate over 
the sea-salt particles found in the MFT . 

According to 3lanchard and Woodcock (1980) , sea-salt 
particles are created during the evaporation of sea water 
drops, and can grow to larger particle sizes due to absorp- 
tion of water under high relative humidity conditions. The 
vertical distribution of sea-salt particles was measured by 
Blanchard and Woodcock (1980) for various wind speeds. The 
presence of the smaller sea-salt particles were found 
throughout the lower free troposphere (i.e., 0. 5-2.0 km), 
and at altitudes of approximately 0.5 to 0.6 km, a sea-salt 
particle inversion was measured. 

The refractive indices for the sea-salt aerosols are 
based on the measurements of Volz (1972) . They are pre- 
sented as a function of wavelength in Table 4. 

The particle size distribution for the maritime aero- 
sol model is a single mode log normal curve. The log normal 
mode parameters, r^ and or, of the distribution were ran- 
domly varied based on real size distribution measurements 
as explained in section 2.1. The mode parameters defining 
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the particle size distributions used in the model are pre- 
sented in Table 3. A normalization of the log normal dis- 

-3 -1 

tributions to 1.08 cm /<m for a radius of 0.1/*m is shown 
in Fig. 4, thus depicting the natural variation in relative 
number concentration one might expect for maritime particle 
size distribution. 

The incorporation of relative humidity effects on dry 
sea-salt particle radius is based on Hanel and Lehmann's 
(1980) mass ratio measurements for sea-salt particles and 
water. The mass ratio measurements are shown in Table 5, 
for increasing and decreasing relative humidities. Also 
shown in the table is the dry particle density and threshold 
water activity, the water activity needed for a particle to 
be completely covered by water. 


Table 5. Mass ratio measurements for sea-salt and 
water for increasing and decreasing rela- 
tive humidity. The dry particle density/* 
and threshold water activity A . are also 
included (Hanel and Lehmann, 1980) . 


MASS RATIO (H 2 0/NaCl) 


rh 

Increasing 

Decreasing 

20 

0.0016 

0.0019 

30 

0.0169 

0.0166 

40 

0.0380 

0.0392 

50 

0.0676 

0.0884 

60 

0.1390 

0.1740 

65 

0.2900 

1.5500 

70 

0.6600 

1.8100 

75 

2.2000 

2.1900 

80 

2.5800 

2.6800 

W 

n / J ^ 

0.7000 

0.6490 

-3 

™ J V 1 ft *5 



W U 

p (gm/cin ) 1.93 
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Fig, 4. Twenty-five log normal particle ? / 

size distributions for the mari- • / 

time (oceanic) aerosol model. All / f 

size distributions are normalized l 

to 1.08 cm’3 for r=0.1/^m to / 6 

show relative variability of the ! 

particle size distributions. r 




The effects of increasing relative humidity on a 
maritime particle size distribution is shown in Fig. 5 for 
the relative humidities 0, 60, 70, and 80%. As explained 
in section 2.2 , the relative humidity effects on particle 
size is minimal at lower relative humidities. As relative 
humidity increases and reaches the critical relative humidity 
for sea-salt aerosols, the particle size distribution ex- 
hibits a significant shift towards the larger particle sizes. 
For decreasing relative humiuity (from near 100%) , Fig. 6 
demonstrates the deliquescent nature of sea-salt particles 
for the relative humidities 80, 70, 60, and 0%. 

2.4 CONTINENTAL (RURAL) AEROSOL MODEL 

The continental (rural) aerosol differs greatly in 
source, size, and composition from oceanic type aerosols. 
Usually, the tropospheric rural aerosols are considered to 
be background sulfate particles with little or no anthropo- 
genic influences. For this study, the composition of the 
sulfate particles were assumed to be ammonium sulfate 
( [NH^J 2 S 0 ^) , which is the principle constituent of the back- 
ground sulfate species (Volz, 1972; Junge, 1953). 

The refractive indices used for the model were based 
on the measurements of Volz (1972, 1973), and are presented 
as a function of wavelength in Table 4. 

As in the maritime model, a single mode log normal 
curve is used to represent the particle size dis tr ibut ion . 
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Fig. 5. Maritime modeled log normal par- 
ticle size distribution for the 
increasing relative humidities 0, 
60, 70, and 80%. The curves are 
normalized to 4 cm - ^ </ 6/m“^ for 
r=0.1//m. 
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Fig. 6. Maritime modeled log normal 
particle si 2 e distributions 
for the decreasing relative 
humidites 80, 70, 60, and 0%. 
The curves are normalized to 
4 cm" m~l for r=0.yni. 
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The log mode parameters, and cr^, were randomly varied 
based on real measurements of the continental size dis- 
tribution and are presented in Table 2. The range of ex- 
pected natural variation for the continental size distribu- 
tions is shown in Fig. 7, where the log normal distribution 
curves have been normalized to 4.9 cm”^ ^m”^ for a radius of 
0.1 pm. 

Mass ratio measurements of ammonium sulfate and water 
made by Hanel and Lehmann (1980) were used to include the 
effects of relative humidity on aerosol size. The mass 
ratio measurements for increasing and decreasing relative 
humidity over the range 0<f<30 are presented in Table 6. 

Also presented are the dry particle density and the 
threshold water activity. 

The effects of relative humidity on dry continental 
particles are very small for relative humidities less than 
70% (Hanel, 1976). A continental particle size distribu- 
tion at relative humidities 0, 60, 70 and 80% is shown in 
Fig. 8 depicting very little variation in the distribution 
due to relative humidity. For decreasing relative humidity 
(from near 100%) , a small hysteresis effect is evident in 
Fig. 9, where a size distribution for relative humidities 
.80, 70, 60, and 0% is presented. The shift towards a "dry" 
particle size distribution occurs at lower relative humidi- 
ties than that for the maritime aerosol. According to the 

results of Fitzgerald et al. (1982) , sulfate particles k 

(£ 
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Fig. 7. Twenty-five log normal particle 
size distributions for the con- 
tinental (rural) aerosol model. 
All size distributions are 
normalized to 4.9 cm“^ for 
r-0.1/im to show relative 
variability of the particle si 2 e 
distributions. 
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Table 6. Mass ratio measurements for sulfate 

and water for increasing and decreas- 
ing relative humdiity. The dry 
particle density and threshold water 
activity A is also included (Hanel 

and Lehmann, 1980) . 

* 


MASS RATIO (NH 4 ) 2 S0 4 


rh 

Increasing 

Decreasing 

20 

0.0154 

0.0636 

30 

0.0270 

0.0713 

40 

0.0462 

0.2080 

50 

0.0732 

0.3970 

60 

0.1180 

0.6090 

65 

0.1640 

0.9330 

70 

0.2720 

1.0600 

75 

0.6090 

1.2000 

80 

1.0160 

1.4200 


(gm/cm *) 1.87 


0.7470 


0.5760 
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Fig. 8. Continental modeled log normal 
particle size distribution for 
the increasing relative 
humidities 0, 60, 70, and 80%. 
The curves are normalized to 
4 cm“3^um~l for r=0.1^um. 
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Continental modeled log normal 
particle size distribution for 
the decreasing relative 
humidities 80, 70, 60, and 0%. 
The curves are normalized to 


4 cm~3 yu for r=0 „ 1 
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m. 



33 


never return to their original dry size. 



III. AEROSOL OPTICAL PROPERTIES 






An examination of the optical properties of aerosols 
with sizes comparable with the wavelength of the incident 
light is presented. The optical calculations are based on 
the Mie Theory solution for a homogeneous sphere intercept- 
ing a plane light wave. Descriptions of the theoretical 
Mie formulas for particle scattering are well documented in 
the literature (Van de Hulst, 1957; Deirmendjian, 1969; 
Kerker, 1969) , and will not be discussed in this paper. 

Mie scattering calculations are often very expensive 
and require a large amount of computing time; therefore a 
cheaper and faster method of calculating aerosol optical 
properties has been explored. A Mie theory algorithm de- 
veloped by Wiscombe (1979) has been tested and used to in- 
vestigate the backscatter phase function (i.e., phase func- 
tion at 180°), and the extinction efficiency in the small 
size parameter range. Average values of the optical pro- 
perties were calculated over a broad range of refractive 
indices and found to be relatively smooth functions of 
size parameter and refractive index. A spline interpola- 
tion of these Mie computations at key refractive indices 
is used to compute backscatter efficiencies and backscatter 
coefficients rapidly. 
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3.1 MIE MODEL COMPARISON TEST 

The Mie model was tested to determine whether the op- 
tical properties were computed accurately. The first test 
examined the ability of the Mie model to calculate the ex- 
tinction coefficient and backscatter phase function for the 
Deirmendjian haze model M. Identical results were calculated 
by the Mie model within four significant digits for the 
wavelength 0.45^111. 

A comparison between the backscatter calculations 
derived from the figures of Mugnai and Wiscombe (1980) and 
the Mie code was used as a second test. Backscatter effi- 
ciencies at a few randomly picked size parameters for the 
wavelength 0.55/^m and m=1.5 are shown in Table 7. 

Table 7. Comparison between the backscatter 
efficiency results of Mugnai and 
Wiscombe (1980) and Mie code based 
on the Wiscombe algorithm. 


SIZE 

MIE C00E 

RESULTS OF MUGU4I 

PARAMETER 


AN0 WISCOMBE 

2.28 

0.52 

0.52 

5.70 

2.10 

2.10 

9.10 

4.47 

4.50 


Based on the above tests, the Mie code has been found to 
sufficiently calculate aerosol optical properties. 
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3.2 BACKSCATTER PHASE FUNCTION FOR m=1.33±0i 

The backscatter phase function of a homogeneous sphere 
with a refractive index m=1.33±0i was examined. Periodic 
fluctuations at small size parameter intervals have been 
discussed by several authors for the phase function of 
spherical water drops (Bryant and Cox, 1966; Nussenzveig, 
1968; Shipley and Weinman , 1978) . For a polydispersed aero- 
sol, an average of the backscatter phase function over these 
periods can be calculated to give the larger scale structure 
of the phase function. Calculations of the average back- 
scatter phase function at various size parameter increments 
is presented. 

The backscatter phase function has been documented as 
exhibiting a quasi-periodic structure with a period of 
AX^ = 0.83 at size parameters near 500 (Bryant and Cox, 
1966; Shipley and Weinman , 1978) . Fig. 10, taken from 
Shipley (1978) , demonstrates the periodic structure of the 
normalized backscatter phase function for the 500 size 
parameter case. For the case of a water drop, Nussenzveig 
approximated the period AX^ » 0.83 of the backscatter 
phase function by using the equation: 
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Fig. 10. The normalized backscatter phase 
function for in = 1.333 and size 
parameters near 500, showing the 
auasi-periodic structure with 
period Ax, = 0.83. This data 
is an extension of the results 
of Bryant and Cox (1966). 
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where the parameter £2 “ 0.262 for m=1.33±0i. The 
parameter £ 2 has only a small dependence on the size 
parameter. 

The backscatter phase function was calculated for the 

-2 

size parameter range 0<X<100. A resolution of AX=10 was 
used for the computations. This resolution is adequate to 
resolve the details of the backscatter phase function, ac- 
cording to the results of Bryant and Cox (1966). A running 
mean of the backscatter phase function was performed over 
an averaging interval of AX^ = 0.83. Two larger quasi- 
periodic oscillations, AX c = 1.1 and AX^ = 14, were found 
in the backscatter phase function as shown in Fig. 11. 
Similar oscillations were reported by Shipley and Weinman 
(1978) for X = 500. 

Since the Mie calculations are so expensive, various 
size parameter increments were tested to optimize the size 
parameter calculations for low cost and accuracy. Figs. 
12a-d represent the change in the average backscatter phase 
function for different size parameter increments. The most 
efficient increment which preserves the feature of the 
averaged backscatter phase function is shown to be 
AX = 0.03. For example, the peaks of the backscatter 
phase function features at the size parameters x=40, and 
x=63.5 can be seen to decrease for AX= 0.03. 



BACKSCATTER PHASE FUNCTION 



SIZE PARAMETER 

Fig. 11. Running mean of the backscatter phase function 
for - 0.01. The running mean was taken over 
the averaging interval « 0.83. 
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Fig. 12b. 


Running mean of the backscatter phase function 
for AX = 0.03. The running mean was taken over 


the averaging interval AX^ = 


0.83. 






BACKSCATTER PHASE FUNCTION 
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Fig. 12d. Running mean of the backscatter phase function 
for AX = 0.05* The runninq mean was taken over 
the averaging intervalAX^ = 0.83. 
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3. 3 BACKSCATTER COEFFICIENTS FROM SPLINE INTERPOLATIONS 

Mie calculations of the averaged backseatter phase 
functions and extinction efficiencies have been used to 
compute backseatter coefficients. The Mie calculations 
were made for small size parameters C<X<1G0 at expected 
values of refractive indices for continental (sulfate) and 
maritime (sea-salt) aerosol types. The average optical 
properties were established to be a relatively smooth func- 
tion of the expected refractive indices, allowing a spline 
interpolation to be used to compute backseatter efficien- 
cies. Integration of the backseatter efficiencies over the 
particle size distributions for each wavelength were made, 
and thus values for the backseatter coefficients were de- 
termined as a function of wavelength. 

The backseatter phase functions and extinction effi- 
ciencies were averaged using the appropriate averaging 
interval ( AX^) for each expected refractive index. The 
values of AX^ were derived from Nussenzvieg ' s equation 
(Equation 5) and are presented in Table 8. 

The averages of the backseatter phase function and 
extinction efficiencies were computed at 2.5 size parameter 
increments for AX = 0.03. The 2.5 size parameter incre- 
ment was used as the midpoint of the averaging interval 
as shown in the example on the following page. 
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Table 8. Period Ax, of the backscatter 

phase function oscillations for 
various refractive indices. 


m 


1.33 

0.83 

1.35 

0.81 

1.37 

0.7ft 

1.39 

0.76 

1.41 

0.74 

1.43 

0.72 

1.45 

0.70 

1.47 

0.69 

1.49 

0.67 

1.51 

0.66 

1.53 

0.64 

1.65 

0.63 

1.57 

0.62 

1.59 

0.60 
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A contour analysis of the averaged backscatter phase func- 
tion for the expected range of refractive indices is shown 
in Fig. 13. The average backscatter phase function appears 
to be a relatively smooth function of the refractive index 
at these small size parameters. A similar smooth function 
was apparent for the average extinction efficiencies, thus, 
allowing average values of the backscatter phase functions 
and extinction efficiencies to be computed by a spline 
interpolation. 

For a given wavelength and refractive index, a spline 
interpolation of the averaged backscatter phase function 
and extinction efficiency values allows for the computation 
of backscatter efficiencies at small size parameters. Back- 
scatter efficiencies were derived by the expression: 

<V IP * Q e 161 

is the backscatter efficiency, 
is the averaged extinction efficiency, 
is the averaqed backscatter phase function. 

The backscatter coefficients were derived by integrating 
the backscatter efficiencies over the log normal particle 
size distributions as shown on page 48. 


where 
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Fig. 13. Contour analysis of the averaged back- 
scatter phase function at small size 
parameters. The running mean for each 
phase function was taken over the 
averaging interval = 0.83, 0.81, 

0.78, 0.76, 0.74, 0.72, 0.70, 0.69, 
0.67, 0.66, 0.63, 0.62, 0.61, .60 for 

refractive indices of 1.33, 1 .5, ..37 
1.39, 1.41, 1.43, 1.45, 1.47, 1.4-, 
1.51, 1.53, 1.55, 1.59, respectively. 
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where 
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dN 
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is the volume backscatter coefficient, 
is the radius, 

is the particle size distribution. 


( 7 ) 


The backscatter efficiencies derived from the inter- 
polated backscatter phase functions and extinction effi- 
ciencies, were used to compute backscatter coefficients as 
a function of wavelength. The range of wavelengths used 
on the calculations was 0.2<\<-2.0 pm. 
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IV. RESULTS AND DISCUSSION 

Backscatter calculations for two key aerosol classes 
were performed as a function of wavelength. These calcu- 
lations included variations due to the natural variability 
of particle size distribution, refractive index, and rela- 
tive humidity. The backscatter coefficients were calculated 
for refractive indices characteristics of sea-salt and 
sulfate aerosol compositions, each over twenty-five parti- 
cle size distributions and seventeen relative humidities 
varying from 0 to 80%. 

Optical properties of atmospheric aerosols are 
weighted by the aerosol size distribution, and aerosol 
scattering is very sensitive to particle radius size 
variability (Deirmendjian, 1969). Aerosol size distribu- 
tions vary in the atmosphere due to variations in meteor- 
ological conditions, sources, sinks, geographical locations, 
and time of year. For very small aerosol radii (i.e., 
Rayleigh scatterer with r<<l), the scattering pattern is 
symmetric. As particle radii increase, the scattering 
pattern becomes asymmetric with more scattering in the 
forward direction. Once the particle radii exceed the 
incident light wavelength, troughs and ridges begin to 
appear in the scattering pattern and most of the energy 
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is scattered into the forward and backward directions 
(Paltridge et al., 1976). Thus, aerosol volume back- 
scatter is very sensitive to size, and it is therefore 
sensitive to natural variability of the particle size dis- 
tribution. Figures 14a and 14b demonstrate this natural 
variability, showing backscatter coefficients for 25 simu- 
lated dry particle size distributions based on available 
measurements from Chapter II. The backscatter curves were 
normalized at a wavelength 0.9yum to highlight the relative 
wavelength dependence. 

The volume backscatter coefficient is assumed to 
follow a power law dependence with wavelength, such that 



This assumption allows a direct calculation of the exponent 
given the backscatter coefficient at two sufficiently 
separated wavelengths. Backscatter coefficients for both 
sea-salt and sulfate aerosol compositions were calculated 
at the wavelengths 0.3, 0.4, 0.55, 0.7, 0.9, 1.05, 1.35, 
1.5, and 1.8/>m. The power law wavelength dependences 
were computed by linear regression to the logarithm of 
these backscatter calculations over four different wave- 
length intervals. The wavelength intervals are 0.3-0. 5 
/im, 0.55-0.9/Vm, 0.9-1.35/im, and 1 . 35-1 . 8 /m, and they 
will be referred to as wavelength regions 1, 2, 3, and 4, 
respectively. 



BACKSCATTER COEFFICIENT (km-1 sr-1) 



Fig. 14a, Maritime (sea-salt) backscatter 
coefficients as a function of 
wavelength for twenty-five dry 
log normal particle size dis- 
tributions. The backscatter co 
efficients were normalized to 
3.7 E-4 (km^sr**) at r=G.9/am 
to highlight the relative wave- 
length dependence. 



BACKSCATTER COEFFICIENT (km-1 sr-1) 



Fig. 14b. Continental (sulfate) backscatter 
coefficients as a function of 
wavelength for twenty-five dry 
log normal particle size dis- 
tributions. The backscatter 
coefficients were normalized to 
1.6 E-4 (km'-lsr"!) at r=0 . 9><m 
to highlight the relative wave- 
length dependence. 
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The theoretical backscatter wavelength dependences 
for sea-salt and sulfate aerosol compositions can now be 
used to distinguish the two compositions* The mean and 
one standard deviation of the £ variations was computed 
in each of the four wavelength regions. If the mean <f 
values are separated by more than the sum of their stan- 
dard deviations, then the two aerosol classes are assumed 
to be distinguishable. Backscatter discrimination between 
sea-salt and sulfate aerosols with the same relative humid- 
ity history is examined in section 4.1. Backscatter dis- 
crimination for sea-salt and sulfate aerosols with 
different relative humidity histories is examined in 
section 4.2. 

4.1 COMPARISON OF SEA-SALT AND SULFATE AEROSOLS AT 
IDENTICAL RELATIVE HUMIDITIES 

The backscatter wavelength dependence of sea-salt 
and sulfate aerosol particles at the same relative humid- 
ity has been compared to determine which wavelength regions 
are best suited to distinguish between the two aerosol com- 
positions. The comparison was made for increasing rela- 
tive humidities from 0 to 80%, with initially dry condi- 
tions, and for decreasing relative humidities from 80 to 
30%, with initially saturated conditions. As relative 
humidity increases, the optical properties of both 
aerosol compositions approaches the characteristics of 
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water and the us-j of S for aerosol discrimination becomes 
less attractive. It is also noteworthy that the wavelength 
regions where sea-salt and sulfate aerosol compositions can 
be discriminated by backscatter is different for increas- 
ing and decreasing relative humidities due to the hystere- 
sis effect of deliquescent materials. 

For the relative humidity increasing from initially 
dry conditions, the water uptake by sulfates and sea-salt 
particles below 50% relative humidities is very small, and 
the wavelength dependence of backscatter does not vary with 
relative humidity. Figures 15-18 show the range of / varia- 
tions for sea-salt and sulfate aerosol compositions under- 
going increasing and decreasing relative humidity varia- 
tions in the four wavelength regions. Each block depicts 
the cf mean (at the center of the block) and one standard 
deviation about the mean. By the definition of distin- 
guishability presented earlier, these figures demonstrate 
that the two aerosol compositions (under dry conditions) 
can be discriminated by cf in the wavelength regions 1 and 
4. The variability of cf at zero relative humidity is due 
solely to variation in the particle size distributions. 

As relative humidity increases to 50% for both aerosol 
compositions, only a small shift in cf is apparent. For 
50% relative humidity, Figures 15 and 18 show that the cf 
discrimination is still possible in wavelength regions 1 
and 4. However, it should be noted that the <f discrimina- 
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Fig* 15, Comparison of the backscatter wavelength dependence 
(/ ) in wavelength region 1 for sea-salt and sulfate 
particles undergoing relative humidity variations. 
Each block represents the mean value (at the center 
of the block) t one standard deviation. 
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Fig. 16. Comparison of the backscatter wavelength dependence 
( ) in wavelength region 2 for sea-salt and sulfate 
particles undergoing relative humidity variations. 
See Fig. 15, 
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Fig. 17. Comparison of the backscatter wavelength dependence 
(c f ) in wavelength region 3 for sea-salt and sulfate 
particles undergoing relative humidity variations. 
See Fig. 15. 


RELATIVE HUMIDITY 


Fig. 18. Comparison of the backscatter wavelength dependence 
(cf ) in wavelength region 4 for sea-salt and sulfate 
particles undergoing relative humidity variations. 
See Fig. 15. 
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tion in wavelength region 1 is approaching the limits of 
the dist inguishability criterion. 

The critical deliquescent point of the aerosol parti 
cles is approached as the relative humidity approaches 60% 
The critical relative humidity of a deliquescent material 
is that relative humidity where a dry salt changes into a 
salt solution. The distribution of <f for both aerosol 
composition starts to shift in all four wavelength regions 
The increase in the average sea-salt aerosol particle 
radius at high relative humidity is much larger than that 
for sulfate aerosol as demonstrated in Figure 3. Thus, 
most of the <5 shifts occur first for the sea-salt composi- 
tion. This is demonstrated in Figure 15, where the sea- 
salt 6 in wavelength region 1 begins to become more nega- 
tive (steeper slope) , while the sulfate <f does not change 
significantly. At 60% relative humidity, sea-salts and 
sulfate aerosols are distinguishable only in wavelength 
region 4 as seen in Figure 18. 

Once relative humidity increases beyond the critical 
relative humidity of a deliquescent material (near 70% and 
above) , that material will begin to acquire the physical 
properties of water, and the ability to distinguish 
aerosol composition by 6 becomes more difficult. Figures 
15-18 show that at these higher relative humidities the 
distributions, except in the sulfate wavelength region 2, 
are converging towards a common range. As aerosol parti- 
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cles become larger they scatter more uniformly with wave- 
length (Van de Hulst, 1957). The shift in wavelength 
region 2 for si/ifates is towards a more positive slope than 
for the J distribution found for sea-salts. These results 
indicate that the o r discrimination between sea-salt and 
sulfate aerosols at higher relative humidities ( 70%) is 
possible in wavelength region 2. 

As relative humidity decreases from saturated condi- 
tions, deliquescent particles will decrease in size at a 
much slower rate than their original growing rate. This is 
due to a deliquescent material’s nature to supersaturate 
with respect to the material in solution. A relative hu- 
midity decreasing from saturated conditions was used in the 
calculations to insure a broad hysteresis loop for each 
composition. Since decreasing relative humidity affects 
particle radius size and refractive index differently than 
increasing relative humidity, the backscatter wavelength 
dependence will not be the same. For sea-salt particles, 
most of the changes in the distribution of <5 occurs at 
relative humidities between 80% and 60%. Once the parti- 
cles dry out at relative humidities less than 60%, the 6 
distribution changes very little and returns to the 
original dry state distribution. For the sulfate parti- 
cles, the shift in <f at higher relative humidities (80- 
50%) is smaller than the shift found for sea-salt. For 
decreasing relative ’.tumidi t ies below 50%, sulfate parti- 
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cles never lose all of the added water. A comparison of 
the sulfate test cases at 30% relative humidity for in- 
creasing and decreasing relative humidity is shown in 
Figure 19 , demonstrating that sulfate particles do not re- 
turn to their original dry cf distributions. The shifts 
found in the two aerosol <f distributions are characteris- 
tic of their hysteresis curves presented in Figure 3. 

For relative humidities decreasing to 50% , the 
ability to discriminate between the aerosol compositions , 
sea-salts and sulfates, within one standard deviation of 
their cT means is possible only in the wavelength region 2. 
This is quite different from the results for increasing 
relative humidity, where discrimination in wavelength 
region 2 was impossible until 70% relative humidity. The 
wavelength region difference between the increasing and 
decreasing relative humidity test cases is related to the 
extra water that both aerosol compositions contain at de- 
creasing relative humidities. As relative humidity de- 
creases to 40%, a shift occurs in the sulfate cf distribu- 
tion towards a more negative slope in the wavelength 
region 2. The two aerosol compositions are no longer 
distinguishable in any of the four wavelength regions. As 
relative humidity continues to decrease, the sea-salt <f 
values return to their original dry distributions . 

Figure 16 show?; the sulfate S distribution shifting toward 
those of sea-salt in wavelength region 2. The cf distribu- 
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Fig. 19. Comparison of the backscatter 
wavelength dependence cf of 
sulfate aerosols at an increas- 
ing relative humidity 30% and * 

a decreasing relative humidity 
30%. See Fig. 15. 
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tions for both aerosol compositions shift away from each 
* other in wavelength regions 1 and 4, thus allowing dis- 

crimination between sea-salt and sulfate aerosols for 
decreasing relative humidity at 30%. 

4.2 SEA-SALT AND SULFATE AEROSOLS AT DIFFERENT RELATIVE 
HUMIDITIES 

The theoretical backscatter wavelength dependence of 
sea-salt and sulfate aerosol compositions has been com- 
pared at different relative humidities to determine which 
wavelength regions can best discriminate between the two 
aerosol compositions. The 6 comparisons were made for an 
aerosol composition at a constant dry relative humidity and 
the other aerosol composition undergoing increasing and de- 
creasing relative humidity variations. This comparison 
applies to dry aerosol layers in the free troposphere with 
most aerosols in the mixed layer. Such aerosol layers were 
observed in the lower troposphere during the Global Tropo- 
spheric experiment (Browell and Shipley, 1982). 

A comparison between the wavelength regions that can 
best discriminate between sulfate aerosols and sea-salt 

* aerosols at increasing and decreasing relative humidity is 
presented in Figures 20-23. In section 4.1, it was shown 
that for zero relative humidity, the two compositions were 

* distinguishable by cf in the wavelength regions 1 and 4. 

1 • 

* As relative humidity was allowed to increase for the sea- 

» 

' salt aerosols (from 0 to 50%) , very little water uptake 
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Fig, 23. Comparison of the backscatter wavelength dependence 
( <$ ) in wavelength region 4 for dry sulfate and sea- 
salt particles undergoing relative humidity varia- 
tions. See Fig. 15, 
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by the particles was apparent* Figures 20 and 23 demon- 
strate that for sea-salt particles at 50% relative 
humidity, the wavelength regions where & can distinguish 
the aerosol compositions is identical to the zero relative 
humidity comparison test case ( i . e wavelength regions 1 
and 4 ) . 

As relative humidity increases from 50%, the dis- 
tribution of sea-salt 6 converges toward a common range. 

For the 60% sea-salt test case, shown in Figure 20, 
values in the wavelength region 1 have shifted towards 
those of dry sulfate, and thus only wavelength region 4 
can be used to discriminate the two compositions. Once the 
critical relative humidity for sea-salts is reached ( 70%), 
the range of J in wavelengths region 4 increases. Figures 
20-23 shows that for relative humidities 75% or above, the 
two aerosol cf distributions are no longer distinguishable 
within the one standard deviation criteria for any of the 
wavelength regions. 

For decreasing relative humidity from saturated con- 
ditions, the two aerosols compositions cannot be distin- 
guished bye f for any of the four wavelength regions. As 
relative humidity decreases from 70% to 60%, the spread of 
for the sea-salt aerosol in the wavelength 4 decreases, 
thus allowing the dry sulfates to be discriminated from 
the sea-salt particles by 6 in the wavelength region 4 as 
shown in Figure 23. Once the relative humidity decreases 
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below 60% , the sea-salt particles start to dry out and 
discrimination is possible in the wavelength region 1. As 
relative humidity continues to decrease to 30%, thee/ dis- 
crimination can identify the two aerosols compositions in 
the wavelength regions 1 and 4, 

The wavelength regions that best distinguish the dry 
sulfates from sea-salt particles at low relative humidities 
are identical for both of the increasing and decreasing 
relative humidity cases. This is characteristic of the 
hysteresis curve for the sea-salt aerosol since sea-salts 
return to their original dry size at low decreasing rela- 
tive humidities. 

The second situation "considered was the comparison 
of dry sea-salt aerosols with sulfate particles undergoing 
increasing and decreasing relative humidity variations as 
shown in Figures 24-27. As mentioned in section 4.1, the 
water uptake by sulfate particles is very small at low 
relative humidities. Figures 24 and 27 demonstrate that 
the small water uptake has little effect on the </ dis- 
criminations, since the compositions can be distinguished 
in wavelength regions 1 and 4 at 20% and 50% relative 
humidities. 

As relative humidity increases to 60% and above, 
the sulfate 6 in the wavelength region 3 shifts away from 
that for dry sea-salts, as seen in Figure 26, allowing the 
<f discrimination to be possible in the wavelength regions 
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Fig. 24. Comparison of the backscatter wavelength dependence 
( </ ) in wavelength region 1 for dry sea-salt and 
* sulfate particles undergoing relative humidity 

variations. See Fig. 15. 
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Fig* 26. Comparison of the backscatter wavelength dependence 
(</ ) in wavelength region 3 for dry sea-salt and 
sulfate particles undergoing relative humidity 
variations. See Fig. 15. 
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Fig. 27. Comparison of the backscatter wavelength dependence 
{/) in wavelength region 4 for dry sea-salt and 
sulfate particles undergoing relative humidity 
variations. See Fig. 15. 
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1, 3 and 4. Once the critical relative humidity for sul- 
fates is reached ( 70%) , the d values in wavelength region 
4 approach those of the dry sea-salt. Figures 24 and 26 
show that only wavelength regions 1 and 3 ran discriminate 
the aerosols at 70% relative humidity. As relative humid- 
ity continues to increase to 80%, the sulfate cT values in 
wavelength region 1 shift towards those of the sea-salt, 
whereas the sea-salt <f values in wavelength region 2 shift 
away from the sulfate values. Figures 25 and 26 demon- 
strate these wavelength region shifts and show that the 
two aerosols can be differentiated in the wavelength 
regions 2 and 3. 

Sulfate particles at decreasing relative humidities 
(from near 100%) are compared with dry sea-salt aerosols. 
As relative humidity decreases from 75% to 60%, <T in the 
wavelength region 2 can best distinguish the two aerosol 
compositions as shown in Figure 25. Once the relative 
humidity decreases to 50%, the sulfate d values in the 
wavelength region 1 shift away from those of dry sea-salts 
shown in Figure 24. In wavelength region 2, Figure 25 
shows the d values have become more similar. At 30% rela- 
tive humidity the d values in wavelength region 4 begin 
to shift away from each other, thus allowing the <f dis- 
crimination. 

The S distribution in the four wavelength regions 
that discriminate between sea-salts and sulfates are dif- 
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ferent for decreasing relative humidities than those 
regions for increasing relative humidities. This differ- 
ence is due to the hysteresis effect of sulfate aerosols 
since sulfate particles never return to their original dry 
size for decreasing relative humidities. 

4.3 SUMMARY 

Sea-salt and sulfate aerosol compositions have been 
characterized and identified by their theoretical back- 
scatter wavelength dependence. Comparisons of c f for the 
two aerosol compositions were made at four wavelength 
regions for increasing and decreasing relative humidities. 
A summary of the results is shown in Tables 9a, 9b, and 9c. 
This study has shown that the influence of the natural 
variability of particle size distributions and relative 
humidity are significant. The effects of increasing and 
decreasing relative humidity on <f were found to be charac- 
teristic of the aerosol hysteresis curves. 
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Table 9a. The wavelength regions that the mean 
of sea-salt and sulfate aerosols at 
identical relative humidities are 
separated by the standard deviation 
criteria (i.e., <r ss+ <7^) . The wave- 
length regions 1, 2, 3, and 4 are 
defined as 0.3-0.55, 0.55-0.9, 0.9- 
1.35, and 1 . 35-1 . 8,4cm, respectively. 


Relative 

Humidity 


Increasinc 


Decreasing 


0 

1,4 

20 

1,4 

30 

1,4 

40 

1,4 

50 

1,4 

60 

4 

70 

2 

80 

2 




Table 9b. The wavelength region that the mean 

of dry sulfate and sea-salt aerosols of 
various relative humidities are separated 
by the standard deviation criteria (i.e., 
tf 'ss+^f ) * The wavelength regions 1, 2, 

3, and 4 are defined as 0.3-0.55, 0.55- 
0.9, 0.9-1.35, and 1 . 35-1 . iyUm, respec- 
tively. 


Relative 

Humidity 

Increasing 

Decreasing 


0 

1,4 

— 


20 

1,4 

-- 


30 

1,4 

1,4 

/ 

40 

1,4 

1,4 

i 

> t 

50 

1,4 

1,4 

60 

1,4 

4 


65 

4 

4 


70 

4 

4 


75 

— 



80 


-- 
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Table 9c. The wavelength region that the mean 
of dry sea-salt and sulfate aerosols 
of various relative humidities are 
separated by the standard deviation 
criteria (i.e., + 0" s f) • The wave- 

length regions 1, 2, 3, and 4 are 
defined as 0.3-0.55, 0.55-0.9, 0.9- 
1.35, and 1. 35-1.8/tm, respectively. 


Relative 

Humidity 

Increasing 

Decreasing 

0 

1,4 

-- 

20 

1,4 

— 

30 

1,4 

1,4 

40 

1,4 

1 

50 

1,4 

1 

60 

1,3,4 

2 

65 

1,3,4 

2 

70 

1,3 

2 

75 

1,3 

2 

80 

2,3 

— 
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V. LIDAR METHOD 

A multi-wavelength lidar system provides information 
on the total volume backscatter coefficient as a function 
of wavelength. This study has theoretically shown that 
aerosol backscatter wavelength dependence is uniquely re- 
lated to aerosol particle composition, and therefore, that 
multi-wavelength lidar measurements may potentially be used 
to infer aerosol composition. The derivation of aerosol 
optical properties from multi-wavelength lidar data can be 
difficult and inaccurate (Shipley, 1978; Klett, 1980), but 
under special conditions simple asymptotic solutions are 
available. Given simultaneous lidar returns to two wave- 
lengths, this chapter presents a method for the computation 
of an average aerosol backscatter wavelength dependence. 

5.1 THE DERIVATION OF 

The lidar equation is given by: 


(2 





7 


79 


where 

P^ power returned 

R range 

lidar constant 

ff Ray^ Rayleigh scattering coefficient 

\ lidar wavelength 

Aerosol to molecular back- 
scatter ratio 


S(R) = 
x 




(■ 




>) 



one way optical depth 


4 


ip normalized aerosol backscatter phase function 

-x 

(Tit * aerosol scattering coefficient. 

aer 

X 

The lidar system constant can be eliminated by taking the 
ratio of the return signal at range R so that at some 
normalization range Rq , such that: 


Kti/ _ (u $,(«)) 

r(U£ (i+s^.i) 



( 10 ) 
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Extracting volume backscatter coefficients from lidar 
measurements becomes problematic when optical thickness 
cannot be neglected. Klett (1980) discusses an analytical 
solution for the total scattering cross-section, given a 
power law relationship between total backscatter and optical 
thickness. This solution assumes scattering by aerosols 
only, and it is sensitive to any range variation of the 
aerosol backscatter phase function. Analytic solutions 
which account for scattering by aerosols and air molecules 
are sensitive to the magnitude of the aerosol backscatter- 
ing phase function. 

A simple analytic solution to the lidar equation can 
be obtained when lidar signals are acquired in the limit of 
small optical depth, such that t ~o . This limit is ap- 
proximated by optically thin aerosol layers which often 
reside in the lower troposphere just above the mixed layer. 
Such layering was observed above the tropical marine mixed 
layer during the 1981 Global Tropospheric Experiment 
(Browell and Shipley, 1982) . Given the additional assump- 
tion of negligible aerosol scattering at the normalization 
range Rg, the lidar equation may then be reduced to: 

nu 

X 

AT-* 0 






Ak 
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where 



P(R)»R 2 


Assuming that aerosol scattering follows a power law 
dependence with wavelength such that: 




(12a) 


(12b) 



(12c) 


where \ 0 is a second lidar wavelenath, then the scatterincr 
ratio must also follow a power lav.’ dependence. The aerosol 
volume backscatter coefficient is then 

The scattering ratio is then oiven by 


(1 2d ) 
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where 6 =oC+y . Given negligible optical depth and no 
aerosol scattering at the normalization range Rq , the wave- 
length dependence for backscatter can be obtained directly 
from lidar measurements at two wavelengths, such that: 



where c f denotes a lidar-der ived estimate of the true value 

7 . 


5.2 ERROR ANALYSIS 

Systematic errors ini' which arise due to the as- 
sumptions of negligible optical depth and no aerosol back- 
scatter in the "clean" region at range R Q are now dis- 
cussed. The effects of both assumptions on the lidar in- 
ferred value are computed for 7 « 2.0 and £ =0.5 
given lidar wavelengths at 300 nm and 600 nm. The sensi- 
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fcivity of these systematic errors to the aerosol scattering 
ratio is also examined. 

To estimate the deviation of cT (the lidar-derived 
value) from £ (the true value) as aerosol optical depth in- 
creases but with no aerosol backscatter at range Rq , Eq, 14 
is expanded to: • • 






0 + M R >) 


e - 1 


(15) 


where 



backscatter factor 


extinction factor 


Given the assumption of small optical depth such that 
AT X «1 , Eq. 15 can be approximated by: 
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If the scattering ratio is large (i.e. , $ (R)>>0), then 

A 

the logarithmic terms in Eq. 16 can be approximated by: 


l"jl - 2aTg 




(17a) 


and 



such that 


zat(9-i) + i n f 


( 18 ) 


Equation 18 shows that the systematic variation in £ 
due to finite optical depth is independent of 5 (R) for 
large scattering ratios. This independence is demon- 
strated in Fig. 28, where the systematic variation of 
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* 

Fig. 28. The systematic variation of 6 - T due to an 
optical depth of 0.1 as a function of the 
scattering ratio, assuming no aerosol scatter 
inq at the norma 1 iza t ion range R . 
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( o - S ) for & - 0.5, 2.0 is shown as a function of 

5 (R) with At = 0.1 and ar = 0.1, 1.3. For small S (R) , 

* . X 
o is shown to be very sensitive to the wavelength depen- 
dence of optical depth, and $ underestimates <£ . 

The systematic variation of <J from S as the scatter- 
ing ratio at range R fl becomes significant and is given by: 


/-r= i 




where At is assumed to be negligible. Letting 
S (R)= 5 ( r q)+A, the normalized backscatter terms in the 

Equation can be expressed as: 


1 + fS, 


1 + fS.W ~ 1= l+f5,(R.> 


and Equation 19 is transformed to: 
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Equation 20 shows that the systematic variation of 
associated with S^(Rg) has no dependence on S^(R) when 
A r » 0. Fig. 29 shows that 4 overestimates T when signifi 
cant backscatter is encountered at the "clean" normaliza- 
tion range R Q . 

5.3 THIN AEROSOL LAYER 

The NASA-Langley UV-DIAL system (Browell, et al., 
1983) , was operated over the state of Ohio as part of the 
1980 Persistent Elevated Pollution Episode Northeast 
Regional Oxidant Study (PEPE/NEROS) . Lidar measurements 
were obtained at two wavelengths to create horizontal 
cross-sect ions of aerosol spatial distributions under the 
flight path of the lidar instrument. During a flight on 
July 24 , 1980, a thin aerosol layer was encountered above 
the continental mixed layer as shown in Figs. 30a and 30b. 
In these figures, the lidar backscatter data is presented 
in gray scale where darker shading represents more aerosol 
backseat ter ing. 

The thin aerosol layer located at 2700 m (AGL) and 
first encountered at 40° 02 # N, 82° 48*W is shown emanating 
from deeper convection upwind. In situ measurements were 
made in the vicinity of these cross-sections, but positive 
identification of the aerosols in this layer was not made. 
Yet, the layer is associated with cumulus cloud activity 
at this altitude from 39° 02*N, 84° 39'W to approximately 
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Fig. 30a. 


Time (r,DT) — height (meters) lidar aerosol cross-sections. 
Darker shading represents higher aerosol concentrations. 
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Time (EDT) - height (meters) lidar aerosol cross-sections. 
Darker shading represents higher aerosol concentrations. 
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38° 44'N, 84° 59'W, where mixed layer aerosols are being 
released into the free troposphere. Thus, the thin aero- 
sol layer appears to have originated from continental 
sources. A relative humidity of 65% was computed for the 
aerosol layer from in situ temperature and dew point mea- 
surements made in the vicinity over Dayton, Ohio. 

The raw DIAL signal profiles at the wavelengths 300 
and 600 nm are shown in Figs. 31 and 32, respectively. 

The signals are shown after background has been subtracted 
and a range squared attenuation correction has been ap- 
plied. The ordinate represents range (or decreasing alti- 
tude) in increments of 150 m. <f was computed throughout 
the layer from range 1400 m to 1700 m with Rq equal to 
1400 m under the assumptions of no scattering at Rq and 

A 

negligible optical depth. Values for <f averaged 0.61+0.21 
throughout the thin aerosol layer. The backscatter wave- 
length dependence of sulfate and sea-salt aerosols were 
compared for increasing and decreasing relative humidities 
in Chapter IV. As shown in Fig. 15, the true aerosol 
classes are not totally distinguishable within the one 
standard deviation criterion at 60 and 70% relative 
humidity. However, the spread in 6 values calculated 
theoretically for sulfate aerosols is approximately the 
same as that determined empirically in the thin aerosol 
layer. Although this agreement does not prove that the 
aerosol layer is composed of sulfate particles, it does 
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Fig. 31. Raw DIAL signal profile corrected for 
range squared attenuation at wave- 
length 300 NM (BI01) . The signals 
minus signal background are presented 
as a function of range R. 
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indicate that it is not dominated by sea-salt aerosols. 
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VI. CONCLUSIONS 


This thesis examined the wavelength dependence of 
volume backscatter for two tropospheric aerosol classes; 
namely, sea-salts and sulfates. Model calculations have 
shown that these two key aerosol classes can be charac- 
terized and identified by their backscatter wavelength 
dependence. Also, it was found that the natural variabil- 
ity of aerosol physical properties such as particle size 
distributions, refractive index, and relative humidity had 
a considerable effect on backscatter wavelength dependence. 
The model results show that a multi-wavelength lidar sys- 
tem may be able to identify the composition of aerosol 
layers by their backscatter wavelength dependence at visi- 
ble and near IR wavelengths. 

The backscatter wavelength dependence of sea-salt 
and sulfate particles was computed for twenty-five log 
normal particle size distributions at seventeen increasing 
and decreasing relative humidities. The variation in the 
backscatter wavelength dependence exponent, <f , for both 
sea-salt and sulfate aerosol compositions was found to be 
a function of the characteristic hysteresis curves of the 
aerosols. For decreasing relative humidities, it was 
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shown that the <f values for the sulfate aerosol never re- 
turns to those for the original dry distributions; whereas, 
the sea-salt <f values do return to values characteristic 
of the dry state. 

A simple method to retrieve <f from a two wavelength 
lidar system was derived for the special conditions of 
small optical depth and no scattering at the normalization 
range R Q . Such conditions are encountered often in the 
free troposphere. An error analysis of these assumptions 
revealed that is not sensitive to optical depth at large 
values of the scattering ratio S.(R) . This was not the 

A 

case for small scattering ratios however. For the assump- 
tion of no scattering at the normalization range R q , small 
/ values were found to be grossly overestimated as S^R Q ) 
increased from zero. This overestimation decreases in 
importance as <f becomes large. 

Estimates of <f were computed throughout a thin aero- 
sol layer for lidar profiles at two wavelengths. The com- 
position of the layer appeared to be of continental origin 
with a relative humidity of 65%. A mean and first standard 
deviation of <f for the aerosol layer was compared with the 
65% increasing and decreasing relative humidity model case, 
and found to be approximately the same. Although this 
comparison does not prove that the composition of the 
layer is sulfate, it does suggest that very little or no 
sea-salt particles are present. 



This study has shown that a multi-wavelength lidar 
does have the potential to characterize broad classes of 
aerosol compositions. The optical properties of sea-salt 
and sulfate aerosol compositions have been shown to be 
uniquely related to backscatter wavelength dependence 
under many physical situations. The characterization of 
thin aerosol layer composition by a two wavelength lidar 
was explored and appears possible. Problems can arise 
when aerosol layers become optically thick or when suit- 
ably "clean” regions for lidar signal calibration are not 
available. Under these situations, an analytical procedure 
can be used to obtain backscatter wavelength dependence. 

The analytical inversion of lidar signals for information 
on the backscatter or total extinction coefficient profile 
requires additional information on aerosol optical proper- 
ties at some calibration range R q (Klett, 1980) , or an 
average value over a range interval, .R^< P. < R 2 (Shipley, 
1978) . This information can be obtained using the average 
slope method in the mixed layer, or by an independent 
measurement of total optical depth by passive techniques 
(Norten et al, 1980). 
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